Resubmitted to ApJ: February 11, 2002 

Preprint typeset using LATp^X style emulatcapj v. 04/03/99 
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ABSTRACT 

We study an empirical model for a unified description of the power spectra of accreting neutron stars 
and black holes. This description is based on a superposition of multiple Lorentzians and offers the 
advantage that all QPO and noise components are dealt with in the same way, without the need of 
deciding in advance the nature of each component. This approach also allows us to compare frequencies 
of features with high and low coherences in a consistent manner and greatly facilitates comparison of 
power spectra across a wide range of source types and states. We apply the model to six sources, 
the low-luminosity X-ray bursters IE 1724-3045, SLX 1735-269 and GS 1826-24, the high-latitude 
transient XTE J1118+480, the bright system Cir X-l, and the Z source GX 17+2. We find that it 
provides a good description of the observed spectra, without the need for a scale-free (1//) component. 
We update previously reported correlations between characteristic frequencies of timing features in the 
light of this new approach and discuss similarities between different types of systems which may point 
towards similar underlying physics. 

Subject headings: accretion: accretion disks - black hole physics - stars: neutron stars: oscillations - 
X-rays: stars 
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1. INTRODUCTION 

Since the first detailed timing studies of accreting X-ray 
binaries made with the EXOSAT and Ginga satellites, the 
aperiodic- and quasi-periodic variability and spectral prop- 
erties of low-magnetic-field (B <; 10 10 G) neutron stars 
(hereafter simply neutron stars) and of accreting black- 
hole candidates (hereafter BHC) have been considered in 
a rather different fashion, despite their striking similari- 
ties. In the case of neutron stars, the sources were soon 
classified into two separate classes (the "Z" and the "atoll" 
sources) based on their spectral and timing behavior (see 
van der Klis 1995a and references therein). This classifi- 
cation applied, with few exceptions, to all bright systems; 
the fainter neutron stars, although not exhibiting the full 
complement of phenomenology, seemed to fit in as well. In 
the case of the BHCs, things appeared to be more complex 
and from the beginning a classification as precise as in the 
previous case appeared to be problematic (see Tanaka & 
Lewin 1995; van der Klis 1995a). Despite this complex- 
ity, however, some common ground was found between 
the phenomenologies of neutron stars and BHCs (van der 
Klis 1994a, 1994b). In particular, strong similarities were 
pointed out in the properties and dependence on X-ray 
luminosity and spectral shape of the powerful broad band 
noise seen in both neutron stars and BHCs, as well as 
between some of the QPOs. These similarities were in- 
terpreted in terms of particular source states occurring in 
both neutron stars and BHCs depending, presumably, on 
mass accretion rate. 

With the advent of the Rossi X-ray Timing Explorer 



(RXTE), our knowledge of the properties of the aperi- 
odic variability of accreting X-ray binaries took a substan- 
tial step forward. The power-spectral densities of neu- 
tron stars revealed two previously unknown QPOs with 
frequencies in the kHz range which, in addition to broad 
continuum components, led to power spectra containing up 
to five simultaneous quasi-periodic oscillations (QPOs). A 
number of characteristic frequencies could thus be identi- 
fied for each power spectrum. The kHz oscillations were 
soon found to behave in a rather regular way and their phe- 
nomcnological properties were fitted into the tight scheme 
organizing the neutron-star systems (see van der Klis 2000 
for a review). BHCs showed fewer QPOs, superimposed 
on strong continuum components, which provided a more 
limited number of characteristic frequencies. Even though 
some new QPOs were discovered, such as the "hecto-hertz" 
QPOs (Morgan, Remillard & Greiner 1997; Remillard et 
al. 1999a, 1999b; Homan et al. 2001) and some <100 Hz 
QPOs turned out to be clearly similar to those in neu- 
tron stars, BHC phenomenology proved once again more 
difficult to classify. There are even indications that the 
similarities between neutron stars and BHCs disappear in 
the newly accessible frequency range above ~100 Hz (see, 
e.g., Sunyaev & Revnivtsev 2000). 

Recently, two studies have revisited the above ques- 
tion and identified a number of power-spectral compo- 
nents with frequencies that follow the same correlations in 
both neutron stars and BHCs. Wijnands & van der Klis 
(1999a, hereafter WK99) considered the break frequen- 
cies of the broad-band noise and the centroid frequencies 
of the low-frequency QPOs in the 0.02-70 Hz range and 
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found these frequencies to follow a similar correlation in 
neutron stars and BHCs covering 2.5 decades in frequency. 
Psaltis, Belloni, & van dcr Klis (1999, hereafter PBK99) 
systematically studied the numerous QPOs and peaked 
noise components in the 0.1-1200 Hz range and identified 
two features whose characteristic frequencies follow a cor- 
relation covering three orders of magnitude in frequency. 
For bright neutron star sources, these two features are the 
"horizontal branch oscillation" (HBO) and the lower kHz 
QPO, respectively. For BHCs and low- luminosity neutron 
stars, they are the the low-frequency QPO and a broad 
peaked component at frequencies between 1 and 10 Hz. 

These studies suggest that the same types of variability 
occur in both neutron-star and black-hole sources, which 
can severely constrain theoretical models for their vari- 
ability properties. Two major classes of models exist 
to explain the variable-frequency QPOs in neutron stars: 
the beat frequency models (e.g., Alpar & Shaham 1985; 
Strohmayer et al. 1996; Miller, Lamb, & Psaltis 1998), 
in which the spin-frequency of the neutron star provides 
one of the characteristic frequencies, and the relativistic 
precession models (e.g., Stella, Victri, & Morsink 1999; 
Psaltis & Norman 2001), in which all QPO frequencies 
arise from general relativistic frequencies in the accretion 
disk. The latter class of models do not depend explicitly 
on the properties of the compact object and are therefore 
also applicable to the case of the variable-frequency QPOs 
observed in BHCs. At the same time, discoseismic models 
offer a promising explanation for the "hecto-Hertz" QPOs 
in these sources (e.g., Nowak et al. 1997; Perez et al. 1997; 
see Wagoner 1999 for a review). 

Similarly to the various models for the quasi-periodic 
variability properties of accreting compact objects, various 
attempts have been made for understanding the physical 
origin of the broadband aperiodic variability in both types 
of sources, which in fact constitutes most of the signal in 
the power spectra of BHCs. This noise continuum can 
often roughly be described by an 1// noise with a char- 
acteristic low-frequency break and may, therefore, be pro- 
duced by a large variety of physical mechanisms. Different 
models have considered the contribution of time-variable 
inhomogencitics ("blobs") in the accretion disks (e.g., Bao 
and Ostgaard 1995; Takeuchi, Mineshige, & Negoro 1995), 
short-lived magnetic flares in the accretion-disk coronae 
(Poutanen & Fabian 1999), erratic injections of soft pho- 
tons in a Comptonizing medium (Kazanas & Hua 1999), 
MHD instabilities in the accretion disks (Hawley & Krolik 
2001), etc. 

From the early days of X-ray astronomy it was clear 
that the power spectrum of the BHC Cyg X-l is not 
a pure power law (Terrell 1972; Nolan et al. 1981). 
With EXOSAT it became evident that the power spec- 
tra of BHCs (Belloni & Hasinger 1990a) as well as atoll 
sources (Hasinger & van der Klis 1989) are not smooth, 
but show "bumps and wiggles". These have been var- 
iously described by allowing breaks or exponential cut- 
offs in power-law fits, by adding broad Lorentzian com- 
ponents, and combinations of these approaches (see refer- 
ences above and, e.g., Berger and van der Klis 1998). For 
the low-luminosity neutron stars, recent work by Olive et. 
al (1998) and Barret et al. (2000) revealed power spectra 
that are quite similar to those of BHCs and atoll sources in 
the low hard state. All these power spectra are flat-topped, 



become steeper towards higher frequency and show bumps 
and wiggles. An approach to the description of these power 
spectra of neutron stars and BHCs adopted by a number 
of authors (e.g., Miyamoto et al. 1991; Olive et al. 1998; 
Nowak 2000) includes no power-law or broken power-law 
components in the fits (except perhaps a single steep power 
law at the very lowest frequencies), but uses only a sum 
of Lorentzian components, some of which are broad (with 
coherence factor Q<1). These authors demonstrated that 
the power spectra can be well modeled without the need 
for power-law components; the bumps and wiggles cor- 
respond to "shoulders" fitted with individual Lorentzian 
components. 

A limitation of the traditional approaches is that peaked 
components (QPO) and broad-band components (noise) 
are treated a priori as different, and therefore are fitted 
with intrinsically different models. This applies to both 
BHC and NS systems. In this paper, we expand and elab- 
orate on the description of the power spectra of accreting 
compact objects as the sum of a small number of distinct 
Lorentzian components, proposing it as a model to fit all 
timing features, independent of their Q values. Here, we 
apply it to power density spectra of a sample of sources in 
specific states, i.e. the Low/hard state of BHCs and the 
low-luminosity state of X-ray bursters, but we plan to in- 
vestigate its applicability on a wider range of sources and 
states. We seek to offer a unified paradigm for the fitting 
and interpretation of the power-density spectra of both 
accreting neutron stars and black holes. This will help to 
unambiguously identify similar variability components in 
neutron-star and black-hole sources and thereby to con- 
strain models for quasi-periodic and aperiodic variability 
of accreting compact objects. We apply this method to 
six additional sources: the low-luminosity X-ray bursters 
IE 1724-3045, SLX 1735-269, and GS 1826-24 (some of 
these data have already been presented by Olive et al. 
1998 and Barret et al. 2000), the recently discovered high- 
latitude X-ray transient XTE J1118+480, which is most 
probably a BHC (McClintock et al. 2001; Wagner et al. 
2001; for the X-ray data see also Revnivtsev, Sunyaev & 
Borozdin 2000), the bright and peculiar NS system Cir 
X-l (see Shircy et al. 1996), and the high-luminosity Z- 
source GX 17+2 (see Homan et al. 2002). Finally, we 
discuss the current status of the correlations between tim- 
ing features of accreting X-ray binaries in the light of our 
new approach. 

2. A LORENTZIAN DECOMPOSITION OF THE OBSERVED 
POWER SPECTRA 

Following Nowak (2000), we describe the power spectra 
of both neutron-star and black-hole sources by a sum of 
Lorentzian components, given individually by 

r 2 A 1 

where r is the integrated fractional rms (over — oo to +oo) 
of each Lorentzian and A its HWHM. We also display them 
as vP v vs. v (sec Belloni et al. 1997; Nowak 2000), where 
P v is the power density in rms normalization (e.g., van der 
Klis 1995b) at frequency v and quote as their "character- 
istic frequency" the frequency v max , at which they 
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Fig. 1. — Examples of Lorentzian models. Left panel: P v rep- 
resentation; right panel: vP v representation. Solid lines represent 
Lorentzians with vq=0.1 Hz and different Q values: the spectra la- 
beled a, b, c, d and e correspond to Q=50, 2, f , 0.5 and 0.1, respec- 
tively. The dashed line is a zero-centered Lorentzian with A=0.1 
Hz. 

attain their maximum in vP u , i.e., 

v max = y^ 2 + A 2 = iv Q v/l + l/(4Q 2 ) (2) 

(see Belloni et al. 1997), where Q = vo/2A. This is the fre- 
quency around which the component contributes most of 
its power per logarithmic frequency interval and is a mea- 
sure for the highest frequency covered by the Lorentzian. 
Note that the maximum in vP v does not correspond to the 
centroid frequency v$ of the Lorentzian but v max > vq. 
The difference is small for narrow QPOs, but becomes 
large in the case of broad features. 

In Figure 1 we show a Lorentzian with vq— 0.1 Hz and 
different values of the coherence Q. The left panel shows 
the usual P v vs. v representation whereas the right panel 
shows vP v vs. v. The shift in v max as Q decreases is evi- 
dent. As the half- width A increases at constant vq in order 
to produce a broader feature, the value of A progressively 
takes over from vq in determining the characteristic fre- 
quency Vmax- The logarithmic FWHM in vP v of the fea- 
ture asymptotically approaches an upper limit of log 2 ^^| . 

For very low values of Q, the function approaches a zero- 
centered Lorentzian (Q=0), which in vP v peaks at A. For 
comparison, the dashed line in Fig. 1 corresponds to a 
zero-centered Lorentzian with A=0.1 Hz. 

For fitting very broad features, even those with a cen- 
troid frequency that is formally larger than zero, a zero- 
centered Lorentzian model is often adequate. As a matter 
of fact, this is what a fitting procedure does: if a power 
spectral component is narrow (Q 3>1), such a routine con- 
verges to a value of vq corresponding to the centroid of 
the observed peak, setting A equal to its half-width. If, 
on the other hand, the component is wide (Q <1), the rou- 
tine chooses a value of vq approaching and set A to fit 
the characteristic frequency of the feature (see Fig. 1). Of 
course, if the component is too broad, one never obtains a 
good fit (this would be the case if the power spectra were 
pure power laws). We discuss below how to deal with the 
intermediate cases. 

Traditionally, the term QPO has been restricted to fea- 
tures with Q >2 (see van der Klis 1995a), while less co- 
herent components have been classified as 'noise'. Broad 
peaked components are sometimes called 'peaked noise'. 
As noted by van der Klis (1995b) this classification is ar- 
bitrary; there is a continuous spectrum of intermediate 



cases, and indeed, there is evidence that a component can 
sometimes change from being a clear band-limited noise 
component into a clear QPO or vice versa (Di Salvo et al. 
2000). Noise components can represent characteristic vari- 
ability frequencies just as QPOs do. However, there is an 
ambiguity in how characteristic frequencies are usually re- 
ported for narrow (i.e., QPOs) and broad (i.e., noise) vari- 
ability components. For QPOs, two parameters are usu- 
ally quoted, i.e., the centroid frequency v$ and the FWHM 
in P v , but it is vq which is usually called the "characteristic 
frequency" of the component. For noise components, only 
one parameter is usually quoted, for example one describ- 
ing the width of the component, such as the HWHM A, 
or some power-law break frequency iVeafc- (Simulations 
show that fitting zero-centered Lorentzians with broken 
power laws and vice-versa yields similar values for the A 
of the Lorentzian and the inbreak of the broke n po wer-law, 
but with systematic differences, see Section p^2.) So, in 
intermediate cases, the definition of the characteristic fre- 
quency is ambiguous; it depends on whether the compo- 
nent is considered to be a QPO or noise. 

Based on the discussion above, we adopt the follow- 
ing algorithm for describing the power spectra. (1) We 
fit a model that is the sum of n Lorentzians and in- 
volves no power-law or broken power-law components. 
(2) As the characteristic frequency of each component 
we adopt its peak frequency in the vP v spectrum, i.e., 
v max = y i/q + A 2 . As noted, if a component is narrow, 
its characteristic frequency v max corresponds to the central 
frequency vq of the Lorentzian (see eq. [2]), whereas, if it 
is broad, the characteristic frequency approaches A. For 
very low-Q components it makes little difference to either 
the derived characteristic frequency or the fit to the data 
to set vq = 0. 

Of course, in the standard "lifetime broadening" inter- 
pretation of a Lorentzian as the power spectrum of an 
exponentially damped harmonic oscillator, the Lorentzian 
centroid and width correspond to two different physical 
time scales (respectively, the oscillation period and damp- 
ing time), whereas in the corresponding interpretation of 
a zero-centered Lorentzian the width is a measure for the 
exponential decay time of an exponential shot. Our de- 
scription of both these types of components in terms of 
a single characteristic frequency v max is motivated pri- 
marily by empirical considerations and has no immediate 
interpretation within the damped oscillator picture. We 
note, however, that in interpretations where broad power 
spectral components arise by the superposition of varia- 
tions covering a range of frequencies rather than by life- 
time broadening, v max is just a measure for the highest 
frequency represented in that range, or the one which 
contributes the highest power. Identifying, in an accretion 
disk model, frequencies with radii in the disk would lead 
to a link between v max and the inner radius of the disk 
annulus contributing to the component. More insight 
on the connection to a physical model might come from 
the comparison of the fractional rms values of the differ- 
ent components in addition to their characteristic frequen- 
cies. As discussed in Sect. 5, two important observational 
facts are that the relative contribution of the broad com- 
ponents in low-frequency systems is such as to produce an 
approximate power law distribution, and that the high- 
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Table 1 
Log of the observations. 



Date Number of Exposure 



Obs. 


Start 


End power spectra 


time (s) 


XTE .11118+480 


A 


00-04-13 


00-05-04 


243 


31104 


B 


00-05-06 


00-05-15 


136 


17408 


IE 1724-3045 


C 


96-11-04 


96-11-08 


239 


30592 


D 


97-02-25 


97-03-25 


70 


8960 


E 


97-04-18 


97-05-20 


59 


7552 


F 


98-02-11 


98-02-15 


236 


30208 


G 


98-02-19 


98-02-20 


383 


49024 


H 


98-02-26 


98-02-28 


243 


31104 


I 


98-03-16 


98-03-18 


147 


18816 


J 


98-04-08 


98-04-10 


232 


29696 


K 


98-05-30 


98-05-31 


119 


15232 


L 


98-06-04 


98-06-04 


95 


12160 


SLX 1735-269 


M 


97-10-10 


97-10-12 


254 


32512 


GS 1826-24 


N 


97-11-05 


97-11-06 


348 


44544 


O 


98-06-07 


98-06-08 


144 


18432 


P 


98-06-12 


98-06-12 


102 


13056 


Q 


98-06-23 


98-06-23 


54 


6912 


R 


98-06-24 


98-06-25 


87 


11136 


GX 17+2 


S 




sec Homan et al. (2002) 




Cir X-l 


T 


96-03-12 


96-03-12 


32 


2048 


U 


96-03-16 


96-03-16 


19 


1216 



est frequency component in BH systems is considerably 
suppressed with respect to NS systems. 

Our choice of components and parameters poses the 
problem of comparison with previous results. In this pa- 
per, we apply conversions from vq and A into v max only 
for the data from Nowak (2000) . All other data discussed 
when dealing with global correlations will have to be ex- 
amined for complete consistency, but since most of them 
involve kHz QPOs, which have a rather high Q, it follows 
that v max m vo, so the effects are small. The Lorentzian 
components of the two power spectra in Nowak (2000) with 
reasonably constrained parameters have Q values rang- 
ing from 0.6 to 2.0, corresponding to correction for their 
characteristic frequencies of 1.9 to 1.03 respectively (see 
eq. [2]). As will be discussed below, we believe that fitting 
all components with Lorentzians will eventually facilitate 
comparison across different types of sources, as it allows 
us to treat broad noise components and QPOs in formally 
the same way without having to decide in advance how to 
interpret each component. This makes it easier to compare 
components in different sources and allows us to make the 
connection with the fits to power spectra of neutron-star 
sources with much higher characteristic frequencies. 

3. OBSERVATIONS 

In order to test our approach on power spectra con- 
taining complex broad-band components, we analyzed 
RXTE/PCA observations of one black- hole candidate and 
five neutron-star systems, covering a wide range of char- 
acteristic frequencies. The black-hole candidate is the re- 
cently discovered X-ray transient XTE J1118+480 (Remil- 
lard et al. 2000, Revnivtsev, Sunyaev & Borozdin 2000). 
This source was in a low hard state during these observa- 
tions The neutron-star systems are three low-luminosity 



X-ray bursters (see Barret et al. 2000): IE 1724-3045 (see 
also Olive et al. 1998), SLX 1735-269, and GS 1826-24, 
the peculiar system Cir X-l (see Shirey et al. 1996), and 
the bright Z-source GX 17+2 (see Homan et al. 2002). 
The three bursters are usually considered to be atoll 
sources which are always in the low (island) state. The 
observation log is reported in Table 1. 

For each observation, usually consisting of several con- 
tinuous intervals corresponding to separate RXTE orbits, 
we produced a series of power spectra from data segments 
lasting 64 or 128 seconds (depending on the total length 
of the observation) and a time binning of 1/1024 seconds, 
normalized after Leahy et al. (1983). These spectra were 
then averaged together and logarithmically rebinned. The 
contribution of the Poissonian statistics was removed (see 
van der Klis 1995a) and the spectra were converted to 
squared fractional rms (Belloni & Hasinger 1990b; van der 
Klis 1995b). The average power spectra from different ob- 
servations were compared and, if they were compatible, 
they were further averaged together. The total result- 
ing numbers of power spectra averaged can be found in 
Table 1. Notice that our dataset B for XTE 1118+480 
is the same as that analyzed by Revnivtsev, Sunyaev, & 
Borozdin (2000), and that our dataset C for IE 1724-3045 
corresponds to the data presented by Olive et al. (1998). 
Also, the power spectra from datasets M (SLX 1735—269) 
and N (GS 1826—24) have been presented in Barret et 
al. (2000). To each of our final average power spectra 
we fitted a combination of Lorentzians as described in the 
previous section. For the fits we used XSPEC VI 1.0. 19. 
The only exception to the procedure outlined above is the 
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Fig. 2.— Left panel: model decomposition for XTE 1118+480. 
Right panel: the same four Lorentzian components but with charac- 
teristic frequencies corresponding to the high-frequency end of the 
WK99 and PBK99 correlations and narrower widths for the two 
high-frequency components. 

dataset for GX 17+2, for which we used one of the power 
density spectra from Homan et al. (2002) for the fits. 

In practice, it turned out that all power spectra (with 
the exception of that for GX 17+2, see below) could be fit- 
ted with the following components (see Fig. 2, left panel): 

• A zero-centered low-frequency Lorentzian L& fitting 
the low-frequency end (the flat top) of the band- 
limited noise visible in all power spectra, with char- 
acteristic ("break") frequency v b . 

• Two zero-centered Lorentzians and L u fitting the 
high-frequency end of the band-limited noise, with 
("lower" and "upper") characteristic frequencies 
and v u . 

• One or two Lorentzians fitting the region around 
the low-frequency QPO. The profile of this QPO is 
sometimes more complex than a simple Lorentzian, 
consisting of a relatively narrow core hpp (the 
"LF QPO") at characteristic frequency vlf, and a 
broader "hump" component L^ with characteristic 
frequency V] x . 

It is evident that, in this description, the band-limited 
noise is fitted with three broad (and in our case al- 
ways zero-centered) Lorentzians. For the lowest-frequency 
Lorentzian, a flat top is actually seen in the power spectra, 
while the true behavior at low frequencies of the other two 
components is masked by components that are more pow- 
erful there. The lowest-frequency Lorentzian, Lb, produces 
the break at Vb in the band-limited noise, the highest- 
frequency one, L u , produces the second break at v Ul and 
the middle one, Lg, essentially "fills the hole" between the 
two, so that between v\, and v u a relatively smooth, ap- 
proximately 1//, spectrum is created, with some "bumps 
and wiggles", as described by, e.g., Hasingcr and van der 
Klis (1989) and Bclloni and Hasingcr (1990a), and on 
which the LF QPO is superimposed. We found that al- 
ways v\) < vlf < Uh < vi < i/ u ; when both the hump 
L/i and the LF QPO Lpp are present, their centroid fre- 
quencies v\ and v\ F typically coincide, but as the hump 
is broader, usually vpp < v^. The three components to- 
gether describing the noise, Lb, L^ and L„, were always 
required by the data, although v u was not in all cases 
well-constrained. Of the two components used to describe 
the low-frequency QPO region, Llf and L/j, at least one 



was always required. Notice that the frequencies for Llf 
and L/j presented in Tables 2 and 3 are centroid frequen- 
cies and HWHM, which are the values obtained by fit- 
ting the data: for these two components the characteristic 
frequencies discussed in the text have subsequently been 
computed by applying Equation 2. 

For the case of GX 17+2, which has much higher char- 
acteristic frequencies, the main differences with respect to 
the scheme described above are: (a) the Lj, component 
is very broad but not zero-centered; (b) the ~Llf compo- 
nent shows two additional peaks corresponding to a sub- 
harmonic and a first overtone (see Homan et al. 2002); 
(c) the two higher-frequency components L^ and L u are 
in this case narrow kHz QPO peaks, as expected; (d) two 
more broad components are needed at frequencies below 



3.1. XTE J1118+480 

The power spectra of XTE J1118+480 was fitted with 
four components (see Fig. 3), the LF QPO already re- 
ported by Revnivtsev, Sunyaev & Borozdin (2000), and 
the three zero-centered Lorentzians Lf,, L^, and h u for the 
band-limited noise. The best-fit parameters are reported 
in Table 2. All components are significant. The values of 
X 2 are high, but an examination of the residuals shows no 
systematic trend. The Q value for the LF QPO is around 
7. 

XTE Jl 118+480 




0.01 0.1 1 10 100 



Frequency (Hz) 

Fig. 3.— Power spectra in vP v form for XTE J1118+480. Lines 
mark the best-fit model and its components. 
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Table 2 

Best-fit characteristic frequencies for XTE 1118+480. For L lf centroid frequencies are quoted. 







L 6 


L 


LF 


L e 


L„ 




Obs. 


v b (Hz) 


<f (Hz) 


A (Hz) 


vi (Hz) 


v u (Hz) 


X^d.o.f.) 


A 

B 


0.062±0.004 
0.052±0.002 


0.115±0.003 
0.081±0.002 


0.009±0.004 
0.006±0.003 


1.27±0.04 
1.10±0.03 


25.88±3.61 
27.88±3.09 


97.4(73) 
156.9(73) 



IE 1724-3045 

0.1 ^ — — 

: Obs. H 




0.01 0.1 1 10 100 



Frequency (Hz) 

Fig. 4.— Power spectra in vP„ form for IE 1724-3045 (Obs. H & 
D). Lines mark the best-fit model and its components. 



3.2. IE 1724-3045 

To fit the ten power spectra of IE 1724-3045, five 
Lorcntzian components are necessary. In the first three 
power spectra (C,D,E), the component in the LF QPO 
region is rather broad, with a Q of 0.7-0.8 (the same com- 
ponent can be seen in Olive et al. 1998). In the remaining 
power spectra, an additional narrow peak appears, with 
Q >1, sometimes as high as Q ~ 15. By simple com- 
parison of the spectra, we identify the narrow component 
as the LF QPO Llf, and the broad component as the 
"hump" L/j. In Figure 4 two representative power spectra 
are shown: in the power spectrum of observation H the 
narrow LF QPO component ~Llf is clearly visible while in 
D it is absent. Note that in the vP v representation the 
hump peaks at a higher frequency than the corresponding 
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Fig. 5.— Power spectrum in vP v form for SLX 1735-269 (Obs. 
M). Lines mark the best-fit model and its components. 

narrow LF QPO. The centroid frequencies of these two 
components are the same, so in a P v plot their peaks coin- 
cide. For observation C the narrow LF QPO is marginally 
present, but not required by the fit, while for observation 
E, like D, it is not detected (see Table 3). 

3.3. SLX 1735- 269 

For SLX 1735-269, we produced only one power spec- 
trum, using the same data as Barret et al. (2000). There 
are additional datasets in the RXTE public archive, but 
their total exposure time is not sufficient for the generation 
of good power spectra. The power spectrum is shown in 
Figure 5. A four-Lorentzian model was used to fit the data. 
The best fit parameters are reported in Table 3. As noted 
by Barret et al. (2000), they are remarkably similar to 
those of some observations of IE 1724—3045. Since there 
is no narrow component, we identify the second Lorcntzian 
as the "hump" component L^. 

3.4. GS 1826-24 

We obtained five power spectra of GS 1826—24. The 
first one (observation N) corresponds to that presented by 
Barret et al. (2000). The two X-ray bursts present in the 
data of observation N (Barret et al. 2000) as well as four 
more X-ray bursts present in the remaining datasets were 
excluded from the analysis. The "hump" component was 
present in all spectra; in four of the power spectra a nar- 
row LF QPO was also observed. Notice that in this source 
there are significant differences between the centroid fre- 
quencies of Lh and Llf- In observations Q and R only 
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Table 3 

Best-fit characteristic frequencies for IE 1724-3045, SLX 1735-269 and GS 1826-24. For L lf and L h 

CENTROID FREQUENCIES ARE QUOTED. 



Lb 



Lz,f 

"LF (Hz) 



Obs. 



"h (Hz) 



A (Hz) 



(Hz 



IE 1724-3045 



A (Hz) 



vt (Hz) 



(Hz) 



x'(d.o.f.) 



— 

D 
E 
F 
G 
H 
I 
J 
K 
L 



0.195±0.010 
0.626±0.038 
0.331±0.026 
0.150±0.008 
0.165±0.012 
0.124±0.006 
0.170±0.009 
0.180±0.008 
0.221±0.017 
0.193±0.014 



0.87±0.03 
2.91±0.13 
1.59±0.09 
0.62±0.02 
0.72±0.06 
0.51±0.01 
0.70±0.03 
0.79±0.03 
1.03±0.07 
0.85±0.07 



0.64±0.03 
1.81±0.29 
1.15±0.12 
0.57±0.02 
0.62±0.05 
0.49±0.02 
0.61±0.03 
0.66±0.03 
0.87±0.07 
0.79±0.07 



10.64±0.51 
13.50±3.80 
14.99±2.14 
9.91±0.37 
9.70±0.52 
8.05±0.20 
9.46±0.52 
10.72±0.44 
11.66±1.01 
11.77±0.79 



191.20±29.30 
256.60±64.50 
198.70±53.70 
200.30±27.25 
166.70±27.18 
142.40±11.88 
203.20±34.89 
208.50±27.98 
184.10±39.70 
178.50±32.53 



166.0(73) 
62.0(70) 
90.0(70) 
108.6(67) 
97.1(67) 
141.6(67) 
93.9(67) 
119.1(67) 
70.3(67) 
108.0(67) 



0.63±0.03 
0.69±0.03 
0.50±0.03 
0.69±0.02 
0.78±0.02 
1.01±0.02 
0.86±0.02 



<0.06 
<0.19 
<0.05 
<0.02 
<0.33 
<0.10 
<0.30 



SLX 1735-269 



0.87±0.11 0.78±0"TTF 

GS 1826-24 



74.65±33.77 



94.4(73) 



M 



0.179±0.014 



10.40±1.63 



~T3 0.272±0.006 1.13±0.02 <0.07 1.37±0.03 1.07±0.03 

O 0.273±0.009 1.00±0.18 <0.10 1.22±0.04 1.07±0.04 

P 0.441±0.017 1.82±0.03 <2.22 2.05±0.06 1.44±0.06 

Q 0.901±0.056 4.26±0.18 2.13±0.35 

R 0.604±0.038 2.28±0.07 <0.45 1.54±0.15 1.58±0.25 



15.07±0.39 
13.54±0.55 
20.09±1.47 
22.20±4.93 
18.61±2.21 



353.80±46.98 
256.40±46.23 
456.25±82.88 
>593.85 
>662.40 



122.1(70) 
128.9(70) 
72.8(70) 
85.0(73) 
106.7(70) 



GS 1826-24 

0.1 ^ r- 

Obs. N 




0.01 0.1 1 10 100 



Frequency (Hz) 

Fig. 6.— Power spectra in vP v form for GS 1826-24 (Obs. N & 
Q). Lines mark the best-fit model and its components. 

a lower limit could be obtained on the characteristic "up- 
per" frequency v u . Figure 6 shows the two power spectra 
of GS 1826—24 with the highest and lowest characteris- 
tic frequencies. Notice the narrow LF QPO component 
present in the power spectrum of observation N and ab- 



sent from observation Q. The best fit parameters can be 
found in Table 3. 

3.5. CirX-1 

A detailed description of the power density spectra of 
Cir X-l can be found in Shirey et al. (1996, 1998). The 
Llf component in this source was seen to vary from 1 to 
^12 Hz (see also PBK99). We selected two observations 
among those reported by Shirey et al. (1996), one at rela- 
tively low (~1.3 Hz) and one at relatively high (~8.5 Hz) 
characteristic frequencies. In order to avoid effects due 
to the variability of the source, for observation U we se- 
lected only the last RXTE orbit, corresponding to the low- 
est characteristic frequencies. Both power density specta 
can be fitted with three Lorentzian components (see Table 
4 and Fig. 7)): a zero-centered one (L b ), a narrow one 
(Llf), and a broad one at high frequencies (L^). 

3.6. GX 17+2 

In order to include one example of a NS source with 
high characteristic frequencies (and kHz QPOs), we se- 
lected one observation of the Z-source GX 17+2 from the 
sample analyzed by Homan et al. (2002), corresponding to 
S z between 0.5 and 0.6 (middle of the horizontal branch, 
see Homan et al. 2002). Notice that, since the two kHz 
peaks in this observation are quite well separated from 
the other components, our fits below 200 Hz do not affect 
the kHz QPOs, so that for these high-frequency peaks we 
could adopt the values reported in Homan et al. (2002). 
We obtained the best fit with a model consisting of six 
Lorentzian components (see Fig. 8): a broad (Q=0.6) 
component which we identify with Lb, a narrow one with 
a sub-harmonic and a first harmonic (see Fig. 8) identified 
with Llf, and two additional broad (Q < 1) components 
with characteristic frequencies below i/j,, which we indicate 
as L' b and L b '. Notice that the L' b component is not very 
significant (about 3a). The best fit parameters can 
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Table 4 

Best-fit characteristic frequencies for GX 17+2 and Cir X-l. All frequencies are i/ v 

SUB-HARMONIC AND FIRST HARMONIC OF h LF ARE NOT REPORTED HERE, (a): VALUES FROM HOMAN ET AL. 

(2002), NOT FROM OUR FITS. 





Obs. 


u v (Hz) v b n (Hz) v b (Hz) v LF (Hz) 


v t (Hz) 


V% (Hz) x 2 (d.o.f.) 


GX 17+2 


S 


0.87±0.11 3.05±0.29 6.93±1.05 38.10±0.25 


537+23" 


798 + 6" 150.2(138) 


Cir X-l 


T 


0.59±0.03 1.27±0.02 


24.09±3.75 


169.6(151) 


U 


4.77±0.16 8.43±0.04 


109.22±10.09 


153.5(168) 



Cir X-l 




0.1 1 10 100 



Frequency (Hz) 

Fig. 7. — Power spectrua in vP v form for Cir X-l (Obs. T & U). 
Lines mark the best-fit model and its components. 

be found in Table 4. 

3.7. Summary of the fits and comparison with power-law 

models 

The three zero-centered Lorentzians Lb, Lf and L M are 
required in many spectra to fit the band-limited noise. A 
narrow low- frequency QPO ~Llf is detected in the black 
hole candidate, and in some of the neutron star spectra 
(preferentially those with the lower characteristic frequen- 
cies - note that the frequencies are low in the BHC as well) . 
The broad hump component is required in most neutron- 
star spectra, but not in the BHC. In IE 1724-3045, when 
Llf and Lh are simultaneously present, their centroid 



GX 17 + 2 




0.1 1 10 100 



Frequency (Hz) 

Fig. 8. — Power spectrum in vP v form for GX 17+2 (Obs. S). 
Lines mark the best-fit model and its components. Notice that, 
unlike all other cases, this power spectrum is not rms-normalized. 

frequencies are identical, but in GS 1826—24 these fre- 
quencies are usually different . 

In order to compare the quality of our current fits with 
that of other models for broad components, we also fit- 
ted the power spectra first with a combination of a broken 
power law and one Lorentzian for the narrow LF QPO 
peak. This model has 7 free parameters (compared with 
the 9 of our model). We applied this model to two repre- 
sentative power spectra (observations C of IE 1724—3045 
and A of XTE Jl 118+480). The best fit X 2 values were 
1220 (/75 dof) and 478 (/75 dof) respectively, to be com- 
pared with our values 166(/73 dof) and 156(/73 dof). 
Clearly, the multi-Lorentzian model fits much better. A 
model consisting of a twice-broken power law (Belloni and 
Hasinger 1990a) plus a LF QPO, with 9 free parameters, 
gives similar x 2 values as our multi-Lorentzian fit. As 
noted in the advantage of our current approach is that 
it facilitates comparison across source types. 

4. GLOBAL CORRELATIONS 

As PBK99 and Nowak (2000) noted, the correlations ob- 
served between the frequencies of the LF QPO and both 
kHz QPOs in luminous neutron stars can be extended to 
lower frequencies by relations between the LF QPO and 
two broad Lorentzians respectively taking over the role of 
the lower kHz QPO (seen in several low-luminosity neu- 
tron stars and in BHCs) and the upper kHz QPO (seen 



9 



in Cyg X-l and GX 339—4). Our results may fit in with 
this. The frequencies i/f,, or vlFi v ti and v u define a 
model that is formally similar to that for luminous neu- 
tron stars. In Figure 2 the left panel, the four frequencies 
arc those of our fit to the BHC XTE 1118+480. In the 
right panel, they represent a power spectrum of a luminous 
neutron-star Z source (e.g., van der Klis et al. 1997); here 
the high-frequency components have been made narrower. 
The band-limited noise Lb and the low-frequency QPO 
Llf are both common also at high luminosities (van der 
Klis 1995a). A broad peaked noise component ("bump" 
or "wiggle") similar to our Li component is seen in sev- 
eral higher luminosity sources as well (e.g., Shirey et al. 
1998; PBK99) and may be a low-frequency manifestation 
of the lower kHz QPO (PBK99). Finally, L u was first 
found by Nowak (2000) in Cyg X-l and GX 339-4 and ten- 
tatively identified as a low-frequency version of the upper 
kHz QPO. The presence of such a component was already 
noted previously (but not interpreted in the same way) by 
Olive et al. (1998) for IE 1724-3045 (our observation C; 
see also Yoshida et al. 1993). 

The power spectra can be more complex than the sum 
of four Lorentzians (see, e.g., Ford and van der Klis 1998; 
Di Salvo et al. 2001; van Straaten et al. 2000 and our GX 
17+2 power spectrum). In particular, the low-frequency 
QPO often shows multiple harmonics, the second some- 
times being the strongest one (e.g., Belloni et al. 1997). 
In the following we use the frequencies of the strongest 
QPO peaks in the power spectrum (but see Sect. 4.2). 
We try to identify the four components discussed here in 
all source classes, and investigate if their frequencies fol- 
low or extend the correlations of WK99 and PBK99. If 
these components at very different frequencies arc thus 
phenomenologically related, perhaps they are physically 
related. In addition to our results, we include additional 
new points on Cyg X-l and GX 339-4 (Nowak 2000), 
XTE J1550-564 (Homan et al. 2001), 4U 1915-05 (Boirin 
et al. 2000), and 4U 1728-34 (Di Salvo et al. 2001). 



4.1. Correlation between the flat-top power level and the 
break frequency 

In Cyg X-l the break frequency Vb and the power level 
of the flat part of the band-limited noise Pfl a t are anti- 
correlated (Belloni & Hasingcr 1990). From their figure, it 
is evident that P nat oc 1/fb- The high-frequency end of the 
power spectrum remains approximately unchanged. This 
is seen in other BHCs (Yoshida et al. 1993; van der Klis 
1994a; Mcndcz & van der Klis 1997), in the neutron star 
4U 0614+09 (van Straaten et al. 2000), and also clearly in 
our two sources showing significant changes in their char- 
acteristic frequencies (Fig. 9 and Fig. 10). 

The Lorcntzian description of the band-limited noise 
suggests a clue to this effect. From Equation (1), for a 
zero-centered Lorentzian, Pfj a t = r 2 /7rA, i.e., when the 
fractional rms r of the Lorentzian is constant, P nat is pro- 
portional to 1/A, and hence, as v a = 0, to \jv max . Hence, 
the observed effect can be reproduced by allowing the char- 
acteristic frequency of Lb to vary, while keeping its total 
fractional rms constant. The dotted lines in Fig. 10 repre- 
sent the predicted P nat vs. Vb relation for constant r. The 
definition of r used here (Eq. 1) involves integration from 



0.1 



0.01 



10" 



10" 



10 " 
0.1 



0.01 



10 ° el 



IE 1724-3045 




GS 1826-24 



o.oi 



0.1 




Frequency (Hz) 



Fig. 9.— Power spectra of IE 1724-3045 (top panel) and GS 1826- 
24 (bottom panel) from this work. They are overlaid for comparison 
of their relative fractional rms and characteristic frequencies. 

— oo to +oo, so the true values of the rms are a factor y/2 
less. Most points approximately follow such correlations, 
with a difference between sources in normalization (i.e., 
different r). 

So, the total rms of Lb is rather constant in each of our 
sources but different between sources. The neutron stars 
all cluster near an rms of ~20%, whereas the two BHCs 
have values of ^30% and ^40%, respectively, consistent 
with what would be expected if rms correlates with a ba- 
sic compact-object property such as mass. The lower rms 
in neutron-star sources could also result from the contribu- 
tion of a non-varying component from the stellar surface. 



4.2. Correlation between the low-frequency QPO and 
break frequency 

This correlation (WK99) contains data from both BHCs 
and neutron stars, to which we add our data except SLX 
1735-269 (which does not have a LF QPO, see Table 3), 
the two accurate points from Table 1 in Nowak (2000), 
GX 339-4 (full dataset) and Cyg X-l (4-15 keV data), cor- 
rected for the low-Q shift (see Sect. 2). and Q-corrected 
GX 339-4 points from Nowak, Wilms & Dove (2002), for 
which we plot their Q2 component as our 1/j. For Cyg X-l 
we used the first "harmonic" of the peak, which is by 
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Fig. 10. — Correlation between the break frequency and flat- 
top level of the corresponding zero-centered Lorentzian. Filled cir- 
cles are from Belloni & Hasinger (1990), other symbols arc from 
this work. The dotted lines represent lines of constant total rms for 
a Lorentzian model. Numbers along the right-hand axis are r 2 as 
defined in Equation 1, i.e., a factor of y/2 larger than the true noise 
amplitude. 

far the strongest (see Nowak 2000 and Sect. 4). We also 
included data for 4U 1728-34 (Di Salvo et al. 2001), but 
not for 4U 1915-05 (Boirin et al. 2000) or XTE J1550-564 
(Homan et al. 2001), as these authors report parameters 
for the QPOs only. 

The result is shown in Figure 11. The two frequencies 
follow the same trend for over three orders of magnitude, 
but the individual correlations in different sources are sig- 
nificantly different. Various systematic effects affect the 
data. The characteristic frequency of the break depends 
on energy and hence on source spectrum and energy band 
(Belloni et al. 1997). WK99 used broken power laws, while 
we used zero-centered Lorentzians. By simulations, we 
find that this causes over- or underestimates of the break 
frequency by <J 15%, unlikely to explain the largest of the 
offsets we observe. Some discrepancies are consistent with 
factors of 2 or 4, suggesting that that some of the QPO fre- 
quencies are not fundamental. Indeed, the sub-harmonic 
in GX 17+2 lies on the WK99 correlation. 

Some sources show an additional broad component, L/j, 
at a slightly higher characteristic frequency than Llf (see 
Sect. 3), which might represent the non- Lorentzian wings 
of the narrow LF QPO. Its characteristic frequency Vh is 
well correlated with vlf' the corresponding points are 
added in Figure 11. This correspondence between a nar- 
row QPO and a broader component at a similar ccntroid 
frequency may be related to the cases in which a narrow 
QPO appears near the break of a band- limited noise (e.g., 
Morgan, Remillard, & Greiner 1997 in GRS 1915+105; Di 
Salvo et al. 2001 in 4U 1728-34; Wijnands & van der Klis 
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Fig. 11. — Correlations between characteristic frequencies v^p 
and vi of variability in neutron stars and BHCs. Small points are 
from WK99, large symbols are from this work and from Nowak 
(2000) and Di Salvo et al. (2000). The crosses below the main 
relation around 1 Hz show the correlation between and L^jr. 

1999b and Revnivtsev, Borozdin & Emelyanov 1999 in 
XTE J1806-246; Homan et al. 2001 in XTE J1550-564). 
Notice that this component is clearly present in the GX 
339-4 data from Nowak, Wilms, & Dove (2002): when only 
one component is present, it is L/j, when two (or more) are 
present, they are Llf and (see Fig. 11). 



4.3. Correlation between low-frequency QPO and Lg 

Figure 12 shows the correlation of PBK99, together with 
the Nowak (2000), Boirin et al. (2000), Homan et al. (2001) 
Di Salvo et al. (2001) data, and Nowak, Wilms, & Dove 
(2002) data. Most points agree quite well with the corre- 
lation, except two GS 1826—24 points and the Homan et 
al. (2001) points (see below). 

The Nowak (2000) points would be further above the 
correlation without Q-correction. Most of the GX 339-4 
points by Nowak, Wilms & Dove (2002) lie a factor of 2 
above the main correlation, again an indication that the 
frequency adopted for Llf might not be that of the fun- 
damental. For the QPO frequencies of XTE J1550-564 
(Homan et al. 2001) we used the strongest of the low- 
frequency harmonic peaks for Llf and the 185—285 Hz 
QPOs for Lg. These points do not follow the main cor- 
relation of PBK99, but extend its "second branch", but 
points corresponding to the beginning of the outburst (Cui 
et al. 1999) in the original PBK99 correlation, follow the 
main branch. The points of 4U 1728—34 by Di Salvo et 
al. (2001) and of 4U 1915-05 by Boirin et al. (2000) 
also follow partly the main branch and partly the second 
branch. 

Di Salvo et al. (2001) concluded that the feature plot- 
ted on this second branch vs. vg is in fact, in our current 
terminology, and hence the main and second branch 
in PBK99 together represent the WK99 relation, in 
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Fig. 12. — vlf an d v t °f variability in neutron starsand BHCs. 
Small points are from PBK99, large symbols are from this work and 
from Nowak (2000), Di Salvo et al. (2000), Boirin et al. (2000), 
Barret et al. (2000). 

accordance with their observation that Lb changes from 
band-limited noise into a QPO when its frequency in- 
creases. Indeed, if we plot i/f, vs. i>n, the points cluster 
around the extension of the second branch to lower fre- 
quencies (not shown in Fig. 12). 

4.4. Correlation between the frequencies of components 

Li and L u 

To this correlation (Fig. 13) between upper and lower 
kHz QPO, we added new data for 4U 1728-34 (Di Salvo et 
al. 2001) and 4U 1915-05 (Boirin et al. 2000), as well as our 
L u components and those by Nowak (2000). Nowak's BHC 
points (with Q-correction) and SLX 1735—269 follow the 
low-frequency extension of the correlation, but the BHC 
XTE J1118+480 and the neutron stars IE 1724-3045 and 
GS 1826—24 are above it. Perhaps this discrepancy is re- 
lated to the presence of other components in the power 
spectrum in the few 100 Hz range. A broad Lorentzian 
with an aprozimately constant 100 — 200 Hz frequency oc- 
curs in low- luminosity neutron stars (see, e.g., Ford & van 
der Klis 1998; van Straaten et al. 2000; Di Salvo et al. 
2001) and might also occur in these sources (van Straaten 
et al. 2001). 

5. DISCUSSION 

Applying the model described in Sect. 2, we have shown 
that the power spectra of low-luminosity burst sources, of 
a BHC in the low state, of a Z source, and of Cir X-l can 
be decomposed into the sum of a small number (4 — 5) 
of Lorentzian components. We obtain, for each power 
spectrum, four characteristic frequencies: ulf or u^, 
vg, and v u , which include the broad noise components as 
well as the narrow QPOs identified in previous studies. In 
this manner, the identification of similar power spectral 
components between different sources (as done in PBK99) 
becomes easier and free of potentially subjective interpre- 
tation (see, however, Nowak 2000 and §|| for a discussion 



Fig. 13. — Correlations between characteristic frequencies vt and 
v u of variability in neutron stars and BHCs. Small points are from 
PBK99, large symbols are from this work and from Nowak (2000), 
Di Salvo et al. (2000), Boirin et al. (2000), Barret et al. (2000). 

of some complications related to the wings and harmonic 
structure of the low- frequency QPO). 

The presence in some power spectra of the additional 
Lh component needs to be discussed in some more de- 
tail. This component seems to be closely related to Llf- 
Sometimes only one of these two components is observed; 
sometimes both are present. When they appear simulta- 
neously, their centroid frequencies are approximately the 
same in IE 1724-3045 but not in GS 1826-24. As al- 
ready mentioned, there are cases of power spectra in the 
literature (see, e.g., Morgan, Remillard, & Greiner 1997), 
for which the characteristic frequency of a band-limited 
noise component is rather close to that of a low-frequency 
QPO. Moreover, Di Salvo et al. (2001) actually observed 
the transformation of one into the other, which also in- 
volved 'transitory' power spectra where both a narrow 
and a broad component seem to be present simultaneously, 
with the QPO sitting on the break of the noise spectrum. 
It seems likely that the connection between Vh and vlf 
is a similar case. Based on their frequency coincidence, 
the two components may be related physically as well. 
However, unlike Di Salvo et al. (2001), we observe that, 
in the case discussed here, there is a preference for the 
presence of the narrow LF QPO when the characteristic 
frequencies are low, not when they are high. That the nar- 
row LF QPO component tends to disappear in the spectra 
with the highest characteristic frequencies suggests that 
caution is necessary when identifying this LF QPO with 
those seen at higher luminosities. It is possible that, af- 
ter one narrow Lorentzian disappears, another one comes 
into prominence, as observed in 4U 1728-34 (di Salvo et al. 
2001). A more detailed discussion of the characteristics of 
this component in a number of BHCs will be presented in 
a forthcoming paper. 

We have also studied the correlations between the prop- 
erties of the power-spectral components and extended the 
correlations found earlier for a wide range of other accret- 
ing sources. Our results are consistent with those of earlier 
studies, but also indicate some discrepancies, which may 
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arise from the fact that different sources may follow sim- 
ilar but not identical correlations (see also Psaltis et al. 
1998 where this was demonstrated quantitatively for the 
case of kHz QPOs in bright neutron-star sources). Before 
discussing in detail the nature of these discrepancies, a 
more homogeneous analysis using identical models would 
be necessary for all sources. 

Our description of the power spectra in terms of just 
Lorentzian components characterized by a single charac- 
teristic frequency v max puts noise and QPO phenomena 
on equal footing. This approach is in accordance with re- 
cent findings from the analysis of RXTE data of the atoll 
source 4U 1728-34 by Di Salvo et al. (2001), who showed 
that in some cases a broad flat-top noise component turns 
into a narrow, peaked component, therefore becoming a 
"standard" QPO. The fact that v max of components of 
widely differing Q seems to maintain (at least on some 
occasions) a simple relation to other parameters suggests 
that perhaps the power distribution of these components 
is related to the extent of disk annuli whose inner radius, 
setting Vmax, varies in a simpler fashion than their width 
(§PJ). It might also be related to the presence of a lifetime 
broadening mechanism that operates at a characteristic 
timescale related to the dynamical frequencies in the disk, 
such as viscous dissipation (see, e.g., Psaltis & Norman 
2001). Finally, in a shot-noise model of the broad-band 
variability, large changes in the Q value of a component 
can be produced by a change in the correlation time be- 
tween different shots (see, e.g., Vikhlinin, Churazov & Gil- 
fanov 1994), accounting thus for such changes. However, 
as noted in Sect. 2, life-time broadening mechanisms do 
not naturally involve v max as a basic parameter; rather, 
centroid and width reflect aspects of the physical process 
which are a priori independent. 

Independent of the particular physical mechanism that 
is responsible for the various power-spectral compo- 
nents, our results strongly suggest that the low-frequency 
Lorentzian (or flat-topped power-law component) is not 
different from the other components. Indeed, not only can 
narrow QPOs turn into flat-topped components, as dis- 
cussed above, but also the frequency is correlated to 
the frequencies of all narrow QPOs. The latter property, 
together with the fact that in BHCs the low-frequency 
"noise" accounts for almost the entire power in X rays, 
suggests that this component is generated at the same re- 
gion in the accretion flow as the other QPOs (see, e.g., 
Psaltis & Norman 2001). 

Our analysis, with the help of the improved statistics in 
the RXTE/PCA data as compared to previous, less sen- 
sitive missions suggests that no 1// noise component is 
needed for the interpretation of power spectra in BHCs 
and low- luminosity bursters. However, the superposition 
of Lorentzian components does produce an overall power 
spectrum that in neutron stars resembles an 1// distri- 
bution over a range of frequencies, and in the BHC XTE 
Jl 118+480 is clearly steeper, but also roughly follows a 
power law. The relative total rms strengths of the various 
Lorentzians are such that a relatively smooth power spec- 
trum (in the P v vs. v representation) is produced. This 
property needs to be addressed by theoretical models as 
much as the values for the characteristic frequencies do. 

Within the context of our results, we can also discuss 
the recent suggestion by Sunyaev & Revnivtsev (2000) of 



a possible way of distinguishing between systems hosting 
a black hole and a weakly-magnetic neutron star based on 
their power spectra. According to these authors, neutron- 
star systems can show significant power above ~ 300 Hz, 
while BHCs do not. The relevant effect can be seen 
comparing Figure 3 to Figure 4: the power spectrum of 
IE 1724-3045 contains much more power at high frequen- 
cies than that of XTE 1118+480, as was already noted by 
Sunyaev & Revnivtsev (2000) and Revnivtsev, Sunyaev & 
Borozdin (2000) for these two sources respectively. 

Our power-spectral decomposition suggests that the dif- 
ference in the power spectra of BHCs and neutron stars is 
caused by two effects. The first is that all characteristic 
frequencies are lower in the case of the BHC, probably re- 
flecting the mass dependence of the dynamical time scale 
in the accretion flow, as was also noted by Sunyaev & 
Revnivtsev (2000). The second is more unexpected: in 
the power spectrum of XTE 1118+480 the fractional rms 
of the Lorentzian with the highest frequency is consider- 
ably lower with respect to the other components. This 
effect is not easy to understand. However, the fact that 
the component is significantly detected in this BHC, albeit 
weaker than in neutron stars, does seem to suggest that it 
is not a feature that requires a neutron-star surface for its 
generation; apparently a component with these character- 
istics can also be produced by the accretion disk around a 
black hole. 

In conclusion, we suggest that the underlying physi- 
cal mechanisms producing similar components in all these 
sources may be the same and that the scatter in their dif- 
ferent correlations is due to a combination of differences 
in the analysis method and differences in the fundamen- 
tal system parameters. Although not all timing features 
in black-hole and neutron-star low-mass X-ray binaries are 
included in this representation (e.g., the VLFN, the normal 
branch oscillations in Z and perhaps some atoll sources — 
but see PBK — , the ~ 1 Hz QPO in dipping neutron stars, 
the ~ 100 Hz QPOs in atoll sources — but see §4.4 — , and 
the HS power-law noise as well as the very-low frequency 
QPOs in some BHCs), the ones included here are impor- 
tant and their understanding would be a major step for- 
ward in the study of variability in accreting X-ray sources. 
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